INITIAL DATA AND THE FINAL FATE OF INHOMOGENEOUS DUST COLLAPSE 



X 



I. H. Dwivedi and P. S. Joshi 



■ Tata Institute of Fundamental Research 

Q '. Homi Bhabha Road, Bombay 400 005 

! India 



> 

m 
(N 
O 
(N 

as 

u 



Proofs to be sent to: 
Dr. Pankaj S. Joshi 
Theoretical Astrophysics Group 
: T.I.F.R. Homi Bhabha Road 
Colaba, Bombay 400 005 
India 



1 



Abstract 

We examine here the relevance of the initial state of a collapsing dust cloud to- 
wards determining it's final fate in the course of a continuing gravitational collapse. 
It is shown that given any arbitrary matter distribution M (r) for the cloud at the 
initial epoch, there is always a freedom to choose rest of the initial data, namely 
the initial velocities of the collapsing spherical shells, so that the collapse could 
result either in a black hole or a naked singularity depending on this choice. Thus, 
given the initial density profile, to achieve the desired end state of the gravitational 
collapse one has to give a suitable initial velocity to the cloud. We also character- 
ize here a wide new family of black hole solutions resulting from inhomogeneous 
dust collapse. These configurations obey the usual energy conditions demanding 
the positivity of energy density. 

PACS numbers: 04.20.Dw, 04.70.Bw 
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1. Introduction 



One of the important questions regarding the naked singularities arising from 
gravitational collapse is the issue of their genericity. Even if naked singularities 
occurred for reasonable forms of matter satisfying energy conditions etc, if this 
occurrence could be shown in a suitable sense to be non-generic, arising from an 
initial data set of zero measure only, then such a result would open up a new path for 
a proper formulation and possible proof for the cosmic censorship hypothesis which 
prohibits the occurrence of naked singularities. To make any progress towards such 
a goal, it is necessary to understand the role of initial data in determining the final 
fate of gravitationally collapsing massive objects in terms of either a black hole 
or a naked singularity. The end product of such a collapse- namely the spacetime 
singularity- is either completely covered by an event horizon, thus producing a black 
hole in the spacetime, or in the case otherwise it has possible causal connection with 
the outside universe, which would require further investigations. 

A crucial aspect of the gravitational collapse of a compact object is the choice 
of initial data. After the star has exhausted its nuclear fuel, it is well known that 
if the mass of the remnant cloud was large enough, gravitational collapse would 
ensue. At the time when the gravitational collapse starts, one of the sets of initial 
quantities are density, pressures, etc., which describe the state of matter distribution 
completely. Another set of initial data has to be in the form of the distribution of 
initial velocities of the cloud towards the center. These two sets of initial data 
for the collapse would actually arise as the end product of the process of the star 
exhausting its nuclear fuel. 

Our purpose here is to analyze, from such a perspective, the role of initial 
data (defined in terms of physical variables such as the initial density and veloc- 
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ity distributions) towards determining the final fate of a gravitationally collapsing 
cloud. We consider here the case of pressure-free dust collapse described by the 
Tolman-Bondi-Lemaitre models [1,2], which is a large class of collapse scenarios 
including the homogeneous case [3] . These models have been studied widely for the 
occurrence or otherwise of naked singularity [4-8] (for further details and related 
issues, see e.g. [9]). However, the space of initial data which would cause a naked 
singularity or otherwise, and it's genericity, yet remain to be properly understood. 
The initial data in this case is characterized by two free functions on spacetime 
(obtained through the integration of Einstein equations), namely the mass function 
M(r) describing the matter distribution, and the velocity distribution function Vi{r) 
of the spherical shells of cloud at the onset of collapse. We study here the conditions 
causing the occurrence of a black hole or a naked singularity and it is shown quite 
generally that given any distribution of mass M(r) (which is part of the initial data 
to be specified) , there always exists rest of the initial data in the form of a suitable 
choice of the initial velocities Vi{r) of the cloud towards the center, such that the 
final outcome of the collapse could be either one in the form of a black hole or a 
naked singularity. 

Since the specification of initial state of matter and densities at the onset of 
collapse constitutes very important initial data physically, we investigate here in 
what way the choice of any particular regular distribution of matter would afi'ect 
the evolution and the final fate of the system in terms of either of these outcomes. 
For example, in the consideration by Oppenheimer and Snyder [3], the density 
distribution chosen for the (pressure free) dust cloud at the initial and also at all 
later times is homogeneous (i.e. p = p{t)), with initial spatial curvature being 
constant (i.e. initial velocity distribution Vi{r) oc r). As is well known, the collapse 
of such a configuration results in a black hole where the resulting singularity is 
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fully covered by an event horizon. A generic density profile, however, need not be 
homogeneous and would be typically higher at the origin and decreasing away from 
the center. Similarly, even if the density profile is homogeneous, the initial velocity 
of the spherical shells need not necessarily be proportional to r. It is thus essential 
in realistic considerations that the effects of both the inhomogeneities present in 
the initial matter distribution, and the variations in the velocity profiles, be taken 
into account while investigating the evolution of collapse as we shall discuss below. 

2. Dust Collapse 

Inhomogeneous dust collapse in the comoving coordinates (i.e. fx* = is given 
by the Tolman-Bondi metric 

ds^ = -dt" + -^—^dr^ + R^{de^ + sin'edcl)^) (1) 
1 -\- E 

T'=e5l5l e = e{t,r)=^ (2) 



• [M I , , R^/^Gi-ER/M) 
R=-\j^+E^t-t„ir) = (3) 

Here T*-' and e denote the energy-momentum tensor and total energy density re- 
spectively. The dot and the prime denote partial derivatives with respect to the 
coordinates t and r respectively, and we have taken R{t, r) < as we are considering 
collapse scenarios only. G{y) is a strictly real positive and bound function which 

has the range I > y > —oo, and is given by 

f sin-\/y ^l-y \ i^ 

G(y) = [—^2 ^ J fo^ 1 > 2/ > 

G{y)^l for y = (4) 
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The functions M, to{r) and E are arbitrary functions of r. In the context of 
the cosmic censorship conjecture in a physicaUy reahstic coUapse scenario, at the 
onset of collapse the star should have regular initial data (density, pressure, etc.), 
and the spacetime be singularity free at this initial epoch of time. The curvature 
singularity would develop later during the final stages of the gravitational collapse. 
We therefore require that at the onset of collapse, at time t — U, The spacetime 
is non-singular and initial data in terms of density etc. do not diverge. Using 
the remaining coordinate freedom left in the choice of coordinate r we rescale the 
coordinate r at the initial epoch of time t — U as 

R{ti,r) = r (5) 

This then implies from equation (3) that to{r), and the important quantity R' are 
given by 

r^/^G(-Er/M) 

(7) 



R' = ((ry - P)Y + (6 - (r/ - ^P)YiG{pY)){-p + 1)^) 



where we have put 

R M' E' E 

Y — — ri — r— — p = r— p — pir) — —r — 

The time t = to (r) corresponds to i? = which is the singularity of spacetime where 
the area of the shells of matter at a constant value of the coordinate r vanishes, 
and corresponds to the time when the matter shells meet the physical singularity. 
Thus the singularity occurs at the coordinate value r at the time t = to{r) > ti and 
hence the range of coordinates is given by 

< r < oo, -oo < t < to(?') (9) 
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From equation (3) (i.e. = E + M/R), we note that M(r) is interpreted as 
the mass function for the cloud, and is related to the density of the matter at the 
onset of collapse by 

M' 

— = e{U,r) = p{r) (10) 

The function E{r) is interpreted as the energy of the cloud, and is related to the 
initial radial velocity Vi (r) of the particles t — ti by the relation 

E{r) = R{0, rf - M{r)/r = V^{r) - M/r (11) 

For a free fall collapse with R{ti) = 0, Vi{r) = 0. For a spherically symmetric star 
of radius rg, the complete set of initial data includes the state of matter inside the 
star, i.e. the matter density p(r), and the radial velocities Vi{r) ait — ti. Note that 
the Kretchmann scalar K = R^^^^Rij^i for the Tolman-Bondi-Lemaitre spacetimes 
is given by 

M'2 MM' M2 ^ ^ 

As the initial data, including the density, should be regular and singularity free at 
t = ti, this implies that throughout this initial surface (i.e. for all allowed values of 
r), we have 

M(r) = r^g{r) (13) 

where g{r) is a differentiable function. We assume that the weak energy condition 
(e(^5 > 0) is satisfied, and therefore M' > (note that at t = U, e = M'/r"^), and 
consequently R is a monotone increasing function of r i.e. R' > 0. Since M(0) = 
(in the case otherwise, there will already be a singularity at r = on the initial 
surface t = U as seen from equation (12)), and due to energy conditions M' > 0, 
we have M(r) > for r > 0. 
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As we are interested in the gravitational collapse scenario of a star, the energy- 
density e is taken to have a compact support on an initial spacelike hypersurface, 
and the Tolman-Bondi-Lemaitre spacetime is matched at some r = const. = Vc to 
the exterior Schwarschild metric, 



The value of the Schwarschild radial coordinate is Vs = R{t, Vc) at the boundary of 
the star. We have M(rc) = 2Ms where Mg is the total Schwarschild mass enclosed 
within the dust ball of coordinate radius r = Tc- Without going into further details 
of the matching conditions, we would like to say a few words regarding the apparent 
horizon. The apparent horizon in the interior dust ball lies at R = M{r). From (3) 
and (6) one can see that the corresponding time t = tah{f) is given by 



Since to{r) > tah for all r > 0, and to{0) = ta/i(0) at r = 0, it follows that, only the 
singularity at r = could be naked, and rest of the singular points for the values 

r > are censored. 

Any radial light ray terminating at this singularity at r = in the past could 
go to the future infinity if it reaches the surface of the cloud r = Tc earlier than the 
apparent horizon dX r = r^- In such a case the singularity would be globally naked. 
The shell-focusing central singularity, which we discuss here, appears at a finite time 
t = to(0),r = 0; however, for g{Qi) = it occurs at an infinite value of coordinate 
time, i.e. to(0) = cxd, unless E{r) = r'^fo{r) in such a way that /o(0) > 0. One can 
therefore also choose suitable radial velocity Vi{r) ( for example, £^(r) oc r", a < 2) 
at the initial t = ti, such that for a given initial density distribution the singularity 
does not form at r = (i.e. it occurs at an infinite value of the coordinate time t). 



ds^ = -(1 - 2Ms/rs)dT^ + (1 - 2Ms/rs)-^drl + r^d^^ 



(14) 



t = tah{.r) = to{T)-MG{-E) 



(15) 
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In such a case, since any singularity forming at r > is covered, this gives rise to a 
black hole as the end state of collapse. 

As mentioned earlier, the main aim here is to show that for a given density dis- 
tribution one can choose the initial radial velocities of the spherical shells such that 
gravitational collapse could result either in a black hole or a central naked singu- 
larity. We analyze here mainly the nature of the central shell-focusing singularities 
at i? = 0. However, as seen from (2), there is a density singularity in space-time 
either when i? = 0, or when R' = which corresponds to a shell-crossing singular- 
ity in the spacetime. A relevant question here may be regarding the occurrence of 
such shell-crossing singularities, and their relevance to our conclusions. The point 
is, during the evolution of the collapse, if a shell-crossing singularity occurs before 
the formation of the physical shell-focussing singularity R — (where the shells 
of matter are crushed to a vanishing area), the metric and the coordinate system 
cannot be readily continued beyond such a shell-cross singularity. It may be noted, 
however, that such shell-crosses are generally regarded to be gravitationally weak 
singularities, through which the spacetime may be continued. It is known that in 
the limit of approach to such a singularity, the spacetime curvatures do not diverge 
sufficiently fast in the context of the dust models being considered here [6] , and at 
least a extension of the spacetime exists through such a singularity. It is also 
possible to discuss the extension of spacetime in a distributional sense through such 
a singularity (see e.g. [9] for a further discussion). In any case, for any given density 
distribution (i.e. for given mass function M(r)) we shall consider here only those 
energy functions E{r) (i.e. the initial velocities Vi{r)) for which the dust collapse 
does not give rise to any shell-crossing singularities for r > 0. 
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To make this specific, we have from equations (4) to (8) 



R' = 



yi-pF / Y^G^pY) Y^G^jpY) r]YlG{pY) 

Vy [yr^Wip)^^ G\p) ^+ 2 



(1 



pYG^(pY) 
G{pY) 



Gip) 



) 



where G 



if) 




Y = R/r < 1 and because the functions G{y), G^{y) are strictly increasing positive 
functions, the coefficients of O and rj are positive in the above expression. Hence 
the conditions for no shell-crossings basically turn out to be > for the case 
when r] = {rM' /M) > 0. However, for the dust cloud satisfying the weak energy 
condition and with the mass function M(r) > as pointed out above, we have 
?7 > 0. Thus, to avoid shell-crossing singularities we must have > 0. As seen 
from (8), this amounts to to{r) being a monotone increasing function of r. Thus 
the singularity curve must be an increasing function and successive singular points 
in space form at successively later times. Physically, this means that the shells of 
matter with increasing values of r arrive at the singularity one after the other at 
later and later times without crossing each other in between. In the case of density 
functions which are decreasing as one moves away from the center ( g'{r) < 0), and 
for the class E <0, the condition above reduces to p' >0 [10]. 

3. Causal Structure and Outgoing Trajectories 

For a particular initial data set to develop either into a naked singularity or a 
black hole, one has to analyze the behavior of causal trajectories and categorize dif- 
ferent situations as to when the spacetime permits outgoing nonspacelike geodesies 
from the singularity, or otherwise. We therefore consider the radial null geodesies, 
given by 



dt 
dr 



± 



(17) 
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where the positive sign represents outgoing solutions and negative sign represents 
ingoing ones. We consider here only outgoing geodesies, and recall for completeness 
some notation and terminology from [7], which will be used here. If there are no 
outgoing null trajectories escaping away from the singularity with their past end 
point at the singularity, then the resulting configuration is a black hole in any case. 
Therefore, for the positive sheet of solutions, we have 

dt R' dR I ^dt, , JE + A/X ^^^^^ , , , 

where u = u{r) is a diff'erentiable monotone increasing function of r such that 
■u(r) > for r > 0, u{0) = and, from equation (7) 

^ = (iVu - Pu)X + i&u - iVu - ^Pu)xiG{-PX)){P + ^)^) = H{X, u) 

(19) 

11 Uj 

X = {R/u), rju = —V, Pu = —P 
ru' ru' 



M 

P=iuE/M), A=— , e, = G^ (20) 

U s/uu' 

The function Qu{f) defined above plays an important role in deciding the end 
state of collapse. The singularity at i? = 0, tt = is a singularity of the differential 
equation (15). For a given set of initial data (i.e. M(r), and Vi{r)), if there exists 
a function u = u{r) such that Qu^/P + ^/X has non-zero definite value as r ^ 
along all X = constant directions, then the null geodesies could terminate in the 
past at the singularity, otherwise geodesies do not terminate there. The necessary 
and sufficient condition [7] for the central singularity at i? = 0, r = (occurring at 
a time t = to{0) > U) to have outgoing characteristics terminating in the past at 
the singularity then is that there should exist a real positive value of X — Xq such 
that 

Xo= lim -= lim ^ = U(Xo,0) (21) 

R^O,u^O U R^O,u^Q du 
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If Xq is positive, this implies that geodesies are outgoing with past end point at 
the singularity, while for negative values of Xq they are ingoing. Therefore, for the 
singularity to be naked, a positive real root of the following equation must exist. 



/ JEo + ^ 

V{X) = U{X,0)-X= I I H{X,0)-X = (22) 

where £'0 = E{0), Aq = A(0). If there are no such positive real roots existing, then 
the singularity must necessarily be covered in the sense that there are no outgoing 
null geodesies from it, and the collapse would result into a black hole. On the other 
hand, if such a value X = Xq exists then near the singularity the characteristics 
have the behavior R = Xqu in the (i?, u) plane. In such a case, since Xq > Aq 
and Aq describes the tangent to the apparent horizon R = F at the singularity, the 
singularity is naked at least locally. 

It follows that for a given initial density profile (i.e. for a given M(r)), if one 
could make a suitable choice of the initial velocity function, that is Vi{r), such 
that (22) admits real positive roots, the gravitational collapse would then result 
in at least a locally naked singularity. The behavior of the roots equation (22) is 
mainly determined by the values of initial velocity Vi{r), and its gradient V-{r) at 
r = 0, for a given mass function M(r). Thus, a right choice of the initial radial 
velocity and gradients for the spherical shells in the neighborhood of the center 
r — towards the center could result in a locally naked singularity, or otherwise. In 
other words, for a given density distribution of matter, the velocity of the spherical 
shells in the neighborhood of the center r = (or the energy function E{r)) at the 
onset of collapse determines the final state of the collapse, ending either in a naked 
singularity or a black hole. 

However, for a faraway external observer viewing the collapse, the singularity 
would be visible if and only if there are causal trajectories with their past end point 
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at the singularity, crossing the boundary r = Tc oi the star with dR/du > 0. An 
important question is, given a star of radius r = Tc and total mass M — M{rc), 
would any singular geodesies reach the outside of the star to a faraway observer 
in the Schwarschild region. To be specific, we can ask for a cloud with a given 
value of the mass and size at the onset of collapse, when a locally naked singularity 
will be globally visible. Our considerations above show that for any given initial 
matter distribution for the cloud, it is the local behavior of rest of the initial data, 
namely the corresponding energy function E{r) or the initial particle velocities, 
which decide the final fate of collapse in terms of a black hole, or a locally naked 
singularity. On the other hand, the global visibility of the singularity will depend 
on the global behavior of the corresponding functions. 

To discuss this issue of global visibility, we note that the future behavior of 
outgoing escaping geodesies, which terminate in the past at the singularity, is de- 
termined by the global behavior of free functions in the form of initial radial velocity 
Vi{r) which appear in the geodesic equations. In order that the singularity be glob- 
ally visible for a given density distribution M(r), the global behavior of the initial 
radial velocity of the spherical shells within the cloud needs to be selected at the 
onset of collapse so that at least some geodesies reach the boundary of the cloud 
before the apparent horizon with a positive value of dR/du. For this purpose, we 
consider the paths of these null geodesies. From equation (15), 

dX _ l^dR UiX,u)-X 
du u du u 

The solution of the above gives trajectories of radial null geodesies in the form X = 
X{u). The necessary condition for a null geodesic to terminate at the singularity 
at -R = 0, -u = is that the equation V{X) — must have a real positive root. Let 
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X = Xq he such a simple root of V{X) = 0. By writing the above as 
dX U{X, u)-X {X- Xo){ho -1) + S 



(24) 



du u u 

and integrating, the null trajectories X = X{u) are given by [7], 

X-Xq = Du''°-^ + u^°-^ J Su-^°du (25) 

where V{X) = {X — Xo){ho — 1) + h{X). The function h{X) contains higher order 
terms in X — Xq, i.e. /i(Xo) = {dh / dX) x=Xq — 0. Also, 

'dV{X) 



dX 



(26) 

X=Xo 



and S = S{X,u) — U{X,u) — t/(X, 0) + h{X) is a finite and continuous function 
in the region t < to{r), i.e. for all values {X,u), and D is a constant of integration 
that labels different geodesies. If the singularity is the past end point of these 
geodesies with tangent X = Xq, we must have X Xq as u — > in (25). Note 
that as X ^ Xo,u ^ 0, the last term in equation (25) always vanishes near the 
singularity, regardless of the value of the constant Hq. This is due to the reason 
that as u ^ 0, X ^ Xq, we have S ^ 0). The first term on the right hand side 
of the equation Du^°~^, however, vanishes only if /to > 1- Therefore, the single 
null geodesic described hy D = always terminates in the past at the singularity 
R = 0,u — 0, with X — Xq as tangent. On the other hand, if /lo > 1 a family 
of outgoing singular geodesies terminates at the singularity in the past, with each 
curve being labeled by different values of the constant D. Thus the necessary and 
sufficient condition for the singularity to be at least locally naked is that a real 
positive root of equation (22) exist. 

In the case Hq > 1, the integral curves given by equation (25) do terminate at 
the singularity with a positive definite value of the tangent X — Xq with each null 
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geodesic being characterized by a value of the constant Z), which is determined by 
the boundary conditions at r = rg. For a singular geodesic reaching the boundary 
of the dust cloud u = Uc = u{rc) with X = (Rc/uc) = Xc we have, 

Xc-Xo = Du'^^-^ +u'^°-^ I Su-^'^du (27) 

Juc 

and hence the equation of such a geodesic can be written as 

X-Xo = (Xe-Xo)(— )'^°-'+«'^°-W Su-^°du (28) 

The event horizon is represented by the null geodesic for which Xc = A(rc). For an 
outgoing null geodesic terminating in the past at the singularity, we have X = Xq > 
Aq as the tangent at the singularity, and the curve is ejected into the region R> F 
where dR/du is positive. Therefore, all the geodesies that reach the line r = Vc 
(where the metric (1) is matched with the Schwarschild exterior) with Xc > A(rc) 
(note from equation (18) that dR/du for such geodesies at r = Tc is positive) would 
escape to infinity, while others would become ingoing. Such geodesies from the 
singularity, that reach future infinity, are then given by the equation (28) with 
Xc > Ac- 

An important subcase to note is when the root equation V{X) = has only 
two real positive roots X = X±. In such a case, ho — l>0 either at X = X^ or at 
X — X-. If /iQ — 1 > at X = X- then ho — 1 < at X — X+. In this scenario 
the geodesies would be terminating at the singularity with tangent X = X_, i.e. 
X — > X_, tt — > 0. From equation (23), the behavior of these geodesies near X = X+, 
assuming that p, rj, (3 are finite as r — > oo, is given by 

(X-X+)^(x- (29) 
u 

and thus as X — > X+ the null geodesies attain arbitrarily large values of r, making 
the singularity globally visible. 
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It follows from the above considerations that the occurrence of a black hole, 
or a locally or globally naked singularity, would generally depend on both the local 
as well as global behavior of the functions M(r) and Vi(r). In particular, the 
occurrence of a naked singularity (which may be locally or globally visible), for a 
prechosen distribution of matter, is determined completely by the initial velocity of 
the spherical shells towards the center at the onset of collapse. In the next section 
we examine the evolution of gravitational collapse from such a perspective. 

4. Black Holes and Naked Singularities 

Before examining the general inhomogeneous dust collapse, we first examine 
the special subclass of homogeneous models, where the density distribution at any 
given instant of time has a constant value in space. The Oppenheimer-Snyder 
models [3] of gravitational collapse fall within this class, where the initial density 
p = e{ti, r) = M' /r"^ = Po = const, everywhere in space, that is 

M(r) = Mor^ (30) 

In their original paper, Oppenheimer and Snyder considered the above constant 
density distribution, with the specific choice of the initial velocity for the cloud 
given by, 

Vi{r) = -VoT, Vo = const. (31) 

It is well-known that such a set of initial data leads to a black hole. 

For any given constant density profile, however, the final result of the collapse 
could be different for different choices of initial velocity profiles as we shall discuss 
here. For example, if one chooses the initial velocity differently so that Vi{r) oc 
rVT+or^, a — const.., or in terms of the function E{r) as 

E{r) = -Mor^(eo + eir^ + 70 (r)r^) 
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where eo,ei, are constants and r^7o(r) is a function of r such 70 (0) = 0, then 
such a choice of E{r) leads to a naked singularity in gravitational collapse (see also 
[7]). The constants cq and ei are chosen here such that 



eoMjf VI -eo 2 m^^^ 2 

With such a choice, u = Mgr^ and the root equation (19) becomes 



(32) 



V{X) = 0=>2x'^ + -To{x-l) ^0, X = x'^ (33) 

For the value of Tq chosen as above, this equation has only two real positive roots, 
namely X± = x'^ > 1, such that x+ > X-, which implies X+ > X_. The exact value 
of these roots depend on the exact value of Tq. It follows that ho > 1 for the root 
X = X_ , and hence a family of outgoing null geodesies terminate at the singularity 
with the tangent X = X_. Also, ho < 1 for the other root X = X+. The behaviour 
of the geodesies is given by equation (29). To give a specific example, one can take 
Fo = 28, and then the roots are X+ = 2.443, X_ = 1.512, and /to - 1 = 0.768 along 
the root X = X_ and ho — I = —0.433 along the root X = X+. The boundary of 
the cloud lies at r = Tc (note that Tc >> 2Ms ^ 1 >> Mori). In the neighborhood 
of the singularity the trajectories are described by 

X-X_ = {Xc-X_){ — f°-^ (34) 

Uc 

For the trajectories with Xc < X_, X decreases as u increases, and they move 
towards the apparent horizon which is a straight line given by X = 1 in [R, u) 
plane. For geodesies with Xc > X_, the value of X starts increasing as u increases 
and they move away from the apparent horizon. The geodesic with Xc= I crosses 
the boundary of the cloud at the intersection of this boundary with the apparent 
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horizon. The trajectories with Xc > 1 reach the boundary at an earher time than 
the apparent horizon and the singularity is globally naked. 

We note in the above case of homogeneous dust collapse that the quantities 
and p are given by = r'^G^{p)p', and p = eo + r^(ei + ^o{r)). Therefore, in order 
that shell-crossing singularities do not occur during the collapse, it is sufficient 
to require that p' > 0, or that p be a monotone increasing function of r. This 
implies that while fixing the energy function E{r), the choice of 70 (r") be such that 
3r^(ei + 70 (r)) + r^Yo{r) > for rc > r > 0. 

It follows from the considerations above that within the framework of homo- 
geneous dust collapse itself, there is a wide variety of new collapse solutions which 
result either into a black hole or a naked singularity, depending on the existence 
or otherwise of the real positive roots of the algebraic equation (33). While it is 
well-known that the Oppenheimer-Snyder homogeneous dust collapse with a spe- 
cific choice of velocity profile results into a black hole, we have seen here that a 
wide variety of different choices of particle velocity profiles would also produce the 
Schwarzschild black hole as the end state of gravitational collapse. 

We now consider the future development of general inhomogeneous density 
profiles during the collapse. As pointed out earlier, for the initial data to be regular 
at the onset of collapse the most general mass function has to be of the form 
M(r) = r^g{r), where g{r) is a differentiable function. For physical reasonableness, 
we first consider the cases where the density is a decreasing function of r away from 
the center (hence ^(0) > 0), and as such g'{r) < 0. If 

^3 



lim 



= (35) 



.9{r)-9{0). 

that means g{r) contains terms that are lower then r^. In such a case, we select the 
energy function as that corresponding to the marginally bound case E[r) — 0, and 
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then u is given by 

-1 ^ 

u = 9 ^ -go"" , 90^ g{0) (36) 



V{X) = X- = ^ X = 1 (37) 



We then have 



Hence the singularity is naked and the geodesies terminate at the singularity with 
the tangent X = 1 in the {R,u) plane. Since E = and g'{r) < therefore there 
are no shell-crossing singularities within the cloud of boundary at r = re- On the 
other hand, if equation (35) is not satisfied that means g{r) contains terms of the 
order or higher. In that case, we choose the velocity profile as given by 

E{r) = -eQr^g{r){l + -f{r)r^) (38) 

which implies p{r) = eo(l + r^7(r)), where 7(r) is a differentiable function such 
that ei = 7(0) > and 3 + rY{r) > 0. Then u = r^, and the constants cq and ei 
are chosen such that 

^ 3ei , 1 3G(eo), 15 ^ , , 

Ti = -i^^= ^) > 13 + ^ V3 (39) 

where Mq = g{0). With such a choice the algebraic equation in question reduces to 

V(X) =0^2a;^ + a;^-ri(a;-l) = 0, X ^ Mqx^ (40) 

This has two real positive roots, namely X± = Mqx'^, such that x+ > X- ^ X+ > 
X_ and therefore the singularity is naked. For this class under consideration, since 
g'{r) < 0, we have p' = eor^(3 + rj'ir)) > and the shell-crossings again do not 
occur. Note that while the only 7(0) and cq values determine the nakedness of 
the singularity, it is the global behavior of 7(r) for the range rc > r > which 
determines that no shell- crossings occur. 
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In the above discussion of evolution of inhomogeneous density profiles, moti- 
vated by considerations of physical reasonableness, we assumed decreasing density 
profiles away from center with g'{r) < and g{0) ^ 0. However, we can consider the 
case of a general mass function M(r) = r^g{r) without these assumptions. We can 
then choose the velocity of the collapsing shells as described by the energy function 
E{r) given by 

E{r) = -r''g{r)p{r), G{p) = ^^^{u) (41) 



where u = F{r) = r^g{r), and 71 (tt) is a monotone increasing function of u with 
71 (it) > and 

'^^^^ =Co>0 (42) 



du 



u=0 



We than have 



drf_ 
du 

The equation V{x) — now becomes 



G„ = -^ = 7i,«(«), 7i,u(0) = Co (43) 
du 



V{X) = {l-J-)^^]-X = 



3 

^ 2x'^ + - 3Co{x - 1) = 

This is the same equation as (33), and would have real positive roots if SCq > 13 + 
in which case the singularity would be naked, otherwise it would be covered. 
Note that for the naked singularity to occur, the above inequality is essential and 
it involves only Co which is the first derivative of 71 (w) at tt = 0. The behavior 
of 71 (w) within the cloud Tc > r > is chosen as per the model one wants and 
its global behavior within the cloud is taken such that there are no shell-crossing 
singularities. The dust models would be bounded, marginally bound or unbounded 
(p > 0, = 0, < 0) depending upon the choice of ^/g'yl > |, = |, < | respectively. 
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For nonoccurrence of shell-crossing singularities, first note that energy condi- 
tions imply that rj = rM' /M > for r > and since 71 is monotone increasing 
function, we have 

e = 7i„MVV^>0 (44) 

Therefore, within the cloud for Tc > r > there are no shell-crossing singularities 
during the evolution of gravitational collapse. 

It thus follows that given any density profile for the cloud at the initial epoch, 
a suitable regular choice of the velocity distribution would make the collapse of 
the cloud result into either a black hole or a naked singularity, depending on this 
choice. This characterizes a wide family of new collapse solutions for inhomogeneous 
dust which result into a black hole, apart from the well-known black hole solution 
resulting from the homogeneous dust collapse. 

Although we have given above only particular initial velocity distributions at 
the onset of the collapse for a given density profile, it must be pointed out that 
many more similar velocity profiles would exist which would also ensure the fate of 
the collapse in terms of either a naked singularity or a black hole as desired. We 
have given here a treatment which includes all possible initial density profiles, while 
various specific density profiles have been worked out earlier by several authors (see 
e.g. [9]) for the occurrence of a naked singularity or black hole. 

An interesting choice of a velocity distribution is worth pointing out in the 
context of the discussion here. The null geodesies in the spacetime are given by the 
differential equation (18). The equation involves two free functions, namely M(r) 
and E{r). Thus, for a given M(r) the exact behavior of the trajectories of photons, 
which are the null geodesies) is described by the way E{r) is chosen. Consider for 
example the mass function of the type given above M(r) = r^g{r), and let us choose 
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E{r) which satisfies 

U{a,u) = a, a > 1 

where a is a constant. In general, U defined by equation (18) is an exphcit function 
of E, M and R and an imphcit function of r; choosing R = au and M = u gives the 
above equation. For such an E{r) the nuU geodesic characterized by the equation 
R = aM terminates at the singularity. This geodesic never gets inside the apparent 
horizon and furthermore along the trajectory we have R' > for r > 0, thus 
avoiding the shell-crossing curve also. The geodesic crosses the boundary of the 
cloud avoiding shell-crossings (if there are any), and makes the singularity globally 
naked. 

We will now consider some further details on the issue of global nakedness for 
the naked singularity formation discussed above. In reality, for an external observer 
viewing the collapse, the singularity would be naked if and only if there are causal 
geodesies (with their past end point at the singularity) crossing the boundary of 
the star r = Tc with dR/du > 0, that is they cross the boundary of the star before 
the apparent horizon. In the above discussion, we have given several models for 
an arbitrary density distribution where the singularity is naked. Do every example 
of a locally naked singularity, under all circumstances, lead to global nakedness? 
The answer is certainly in the negative. As such local nakedness already implies 
that for a set of given M(r) and E(r), which leads to a naked singular spacetime, 
causal geodesies which terminate at the singularity in the past do escape into the 
region R > M near the singularity and as such have dR/du positive in a finite 
neighborhood of the singularity, however small. It follows that one can always 
choose the boundary at r = rg > 0, however small, so that in this region and at 
the boundary we have dR/du > 0. Therefore one can always select the total mass 
M(rc), and the size Tc of the cloud, such that the escaping null geodesies would 
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reach Tc with dR/du positive for such examples, making the singularity globally 
naked. Similarly one can, by choosing the boundary suitably, make it globally 
invisible in these cases. Thus for a given set of density and velocity distributions, 
which leads to a locally naked singularity, one can always choose the the set of total 
Schwarschild mass M(rc), and the size of the cloud Tc so as to make the collapse 
either globally naked or otherwise. 

For an arbitrary set of mass M and size rc, let us consider the general mass func- 
tion and the example given by equation (41) where E(r) = —r^g{r)p{r), G{p) = 
y^7i(w), u = M. For the physically reasonable case g{0) ^ one can take the 
bound case p < 0. In this case, between the cloud rc > r > the only con- 
dition one has on the function 71 is that it be a monotone increasing function, 
i.e. 71 = 71, uM' > with 71(0) > 0. At the center at r = one should have 
715U (0) = Co > 13 + ^-\/3, and then the family of geodesies terminate at the sin- 
gularity. The root equation (22) has only two real positive roots > X-. 
The geodesies terminate at the singularity with the tangent X = X_ > 1 in (i?, u) 
plane and attain arbitrarily large value of r at X = In {R,u) plane the ap- 
parent horizon is a straight line R = u and X± > 1 are the roots of equation (22). 
Geodesies with Xc > X_ will escape to the boundary with dR/du positive and 
the singularity would be globally naked. In cases we considered with ^(0) = 0, the 
choice of E{r) corresponds only to the unbound class, however, the conditions on 
71 remain the same for no shell-crossings as well as for the singularity to be naked. 

So far we have concentrated on initial matter distributions of the star. One may 
also ask whether the converse is true, that is, given a regular initial data in the form 
of E{r), whether there always exists a regular M(r) which could produce either one 
of the black hole or a naked singularity. For example, the collapse initiating from rest 
with Viir) — (i.e R — Q initially) falls under such a situation and can be considered 
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similarly by making a suitable choice of the mass function M(r). The answer in 
general seems to be in affirmative, and one can develop the arguments along similar 
lines as above, though we will not go into the details here. Another point to note is 
regarding the curvature strength of the naked singularities considered here. Not all 
the examples considered here need be that of strong curvature type, except the ones 
where the root equation admits two real positive roots, which will be necessarily 
strong. These details will be given elsewhere [11]. 

5. Concluding Remarks 

A relevant question here is whether results such as above could be generalized 
to more general equations of state. In fact, recent work [12] has indicated that the 
phenomena of naked singularity or black hole formation in collapse is closely related 
to the specification of initial data. It has been pointed out there that if a naked 
singular spacetime (or a black hole) exists for a particular equation of state then for 
all sufficiently close equations of state there would also be a naked singular solution 
(black hole) in spherically symmetric collapse. Hence similar conclusions as above 
should hold for such more general equations of state as well. Further details on this 
and related issues will be discussed elsewhere. 

We have shown here that for gravitational collapse of dust clouds with arbitrary 
initial matter distributions, the final fate of collapse in terms of a black hole or a 
naked singularity is fully determined by the choice of rest of the initial data. It 
follows that given any regular matter distribution, its collapse can always lead to 
the formation of a black hole provided the remaining initial data is chosen suitably; 
and that such a choice is always possible. In other words, given the initial density 
distribution, we have given here a specific procedure to fine tune rest of the initial 
data so as to produce a black hole as the end product of gravitational collapse. 
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Such a conclusion no longer needs the assumption of cosmic censorship conjecture. 
Usually, in black hole physics, the truth of cosmic censorship is always assumed in 
the sense that collapse from regular initial conditions is necessarily taken to be a 
black hole. However, no proof or a mathematically suitable formulation of such a 
hypothesis is available so far. We have provided here a more specific characterization 
of initial data for the formation of black holes without such an assumption. Thus, 
these considerations bring out a wide new class of collapse solutions producing a 
black hole from the inhomogeneous collapse of matter. We also pointed out initial 
data configurations resulting into naked singularities for any given distribution of 
matter. In view of such a situation, it appears important to identify and determine 
physically the possible initial data sets responsible for either of these phenomena. 
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